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SUMMARY 


An experimental program was conducted in the Lewis Research Center's 
zero-gravity facility to examine the liquid flow patterns that result from 
the axial-jet mixing of ethanol in 10-cm-diameter containers under zero-, 
reduced-, and normal-gravity conditions. Four tank configurations were used 
for this study: a spherical container; a convex hemispherically ended 

cylindrical container; and two Centaur liquid-hydrogen-tank models, one with 
slosh baffles and one without. Four distinct liquid flow patterns were 
observed: dissipation of the liquid jet in the bulk-liquid region, geyser 

formation, liquid collection at the end of the tank opposite the jet outlet, 
and complete circulation of the liquid along the tank walls. Dimensionless 
parameters were developed that characterized the liquid flow patterns and 
the bulk-liquid mixing phenomena. 


INTRODUCTION 


The use of cryogenic propellants in NASA's space program introduced 
problems of propellant management (ref. I). During missions in space, cryo- 
genic propellant tanks require periodic venting to keep the tank pressure 
within acceptable limits as a result of various heat inputs. Currently, 
orbiting propellant tanlcs - are " vented by using settling rockets to position 
and maint ain the liquid away from the vent and thus to permit the venting of 
only vapor from the tank. 

Efforts to improve the performance of propulsive stages have introduced 
the thermodynamic vent concept, a method of controlling tank pressure other 
than liquid settling and direct venting (ref. 2). Thermodynamic vent system 
concepts are particularly attractive when cryogenic payloads for the space 
shuttle are considered. Under many normal operating conditions and all 
abort modes, payload requirements cannot dictate shuttle operations, so 
settling to relieve the payload tank pressure would be impossible. Thermo- 
dynamic vent systems should impose no operating constraints on the space 
shuttle. 

In the thermodynamic vent system a small amount of the cryogenic propel- 
lant is evaporated to offset the unavoidable heat addition to the propellant 
tank. Cryogenic 1 iquid is withdrawn from the tank and passed through a 
Joul e-Thompson valve with a resultant pressure and temperature reduction. 
This cold, two-phase fluifl is then introduced into a heat exchanger, where 
evaporation continues and heat absorption takes place before the resulting 
vapor is vented overboard. 

For large cryogenic tanks internal thermodynamic vent systems employing 
compact heat exchangers and fluid mixers are preferred because of their 
relatively low weight. However, the internal thermodynamic vent system 
introduces the requirement of circulating or mixing the fluid in the tank so 
that effective fluid cooling occurs and thus reduces the tank pressure. 

Thermodynamic vent systems with internal, compact heat exchangers and 
mixers have been designed, fabricated, and tested in cryogenic fluids under 
normal-gravity conditions (ref. 3). The ability of these devices to control 
tank pressure at 1 g has led to their general acceptance as the planned 
pressure control technique for orbit transfer vehicles and other future 
cryogenic tankage designed for weightless operation. 

Techniques for designing thermodynamic vent systems have not been veri- 
fied experimentally under low-gravity conditions. Specifically, the design 
of the mixer is strongly influenced by the acceleration-dependent fluid 


dynamics necessary to mix the tank contents in a predictable manner (refs. 4 
and 5). The difference between normal gravity and weightlessness will in- 
fluence not only the liquid-vapor configuration before mixing, but also the 
liquid momentum or velocity required to establish the desired liquid flow 
pattern and degree of mixing. 

In earlier studies the author examined the liquid flow patterns that 
resulted when an axial liquid jet (a circular liquid jet directed along the 
axis of the container) was used to mix ethanol in 10-cnvdiameter spherical 
(ref. 4) and cylindrical (ref. 5) containers (fig. 1) under weightlessness. 
Several tests were also conducted under normal-gravity conditions for com- 
parison. The effects of the volume of liquid in the tank, the position of 
the liquid-jet outlet, and the liquid-jet flow rate on the resulting liquid 
flow patterns were examined. Reference 5 presents the author's rationale 




(a) Cylindrical tank with convex (b) Centaur liquid-hydrogen-tank 

hemispherical ends. model without slosh baffles. 

(c) Spherical experiment tank. 

Figure 1. - Experiment tanks. 





for using axial liquid jets to provide fluid mixing and the criteria for 
selecting liquid-jet inflow rates. 

The experimental program was continued by examining the liquid flow pat- 
terns that resulted when an axial liquid jet was used to mix the contents of 
partially filled spherical and cylindrical containers under reduced-gravity 
conditions. Acceleration levels were chosen that simulated the atmospheric 
drag that would exist on a full-scale propellant tank in low Earth orbit. 

Of equal importance with observing liquid flow patterns was examining the 
motion of the ullage, or main vapor region. Vapor encapsulation of the 
simulated pump (i.e., vapor ingested in the region of fluid withdrawal) 
would lead to failure of the liquid circulation concept. 

The objectives of the program were (1) to determine how the acceler- 
ation environment, coupled with the previously studied range of variables, 
affects the resulting liquid flow patterns and (2) to establish analytical 
correlating parameters that would characterize the observed fluid dynamic 
phenomena. 


SYMBOLS 

a acceleration, cm/sec^ 

2 

Bo Bond number, apR .la 

3 

D diameter, cm 

F flow characterization parameter, (X + Y We) / ( 1 + Z Bo) 

G ratio of geyser height to tank radius, h g /R t 

g acceleration due to gravity, 980 cm/sec^ 

h height, cm 

M mixing parameter, Q e/V^ 

Q liquid-jet volumetric inflow rate, cm^/sec 

R radius, cm 

Re Reynolds number, pvD/y 

V volume, cm^ 

v liquid-jet velocity, cm/sec 

2 2 

We Weber number, pv Q R Q loD 

X,Y,Z nondimensional constants 
e bulk-liquid mixing time, sec 

viscosity, g/cm-sec 


u 


p 


density, g/cm 3 
o surface tension, g/sec 3 

Subscripts: 
b bulk liquid 

c completely turbulent 

g geyser 

j liquid jet 

£ laminar 

o outlet 

p partially turbulent 

t tank or container 


APPARATUS AND PROCEDURE 

The Lewis zero-gravity facility was used to obtain the experimental data 
for this investigation. The facility, the experiment package, and the pro- 
cedure for conducting the tests are described completely in the appendix of 
reference 5. 

A schematic drawing of the liquid flow system is shown in figure 2. A 
piston pump, driven by an electric motor and gear mechanism, forces the 
ethanol test liquid out through the central tube and forms an axial liquid 
jet. Simultaneously, liquid is withdrawn from the bottom of the experiment 
tank and fills the region below the piston pumps. Interchangeable 
0.4-cm-diameter liquid-jet outlet tubes of different lengths could be fit 
into the opening at the top of the piston pumps. 

The speed of the piston pump, and thus the liquid-jet volumetric flow 
rate, is controlled by various combinations of gear drives and applied 
voltage to the electric motor. As the piston moves, an electrical contact 
slides along a linear resistor mounted parallel to the piston pump. A 28 -V, 
direct-current circuit, including the linear resistor that acts as a voltage 
divider, provides the input to a telemetry system that transmits to a 
receiver and a strip-chart recorder in the facility control room so that the 
position of the piston pump as a function of time can be determined. 

Because the piston drive mechanism occupies approximately 30 percent of the 
piston sleeve volume, the liquid volumetric withdrawal rate is only 
70 percent of the liquid-jet volumetric flow rate. 

During a test under zero- or reduced-gravity conditions the first 2 sec 
of test time was allocated for the bulk liquid to achieve a reduced- or 
zero-gravity configuration. The liquid flow system was then activated by 
timers for the remaining test time, 3.1 sec. For simplicity, the same 
sequence and timing of operations were used for the normal-gravity tests. 
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lliiilii.l 


.^Jet outlet tube 




— Piston drive mechanism 


Figure 2. - Liquid flow system. 


A high-speed motion picture camera was used to photograph the initial 
bulk-liquid motion and the liquid flow patterns within the experiment tanks 
after the piston pump was activated. The liquid jet was dyed to clarify the 
liquid flow patterns. A digital clock was included in the camera field of 
view so that the elapsed time from the start of the test would be recorded. 

Four 10-cm-d i ameter experiment tank configurations were used for the 
study. The four tanks are identified in table I by the following letters: 

A - convex hemi spherical ly ended cylindrical tank with a length-diameter 
ratio of 2 and a volume of 1310 cm 3 ; B - a Centaur liquid-hydrogen-tank 
model with one convex hemispherical end, one concave ellipsoidal end, no 
slosh baffles, a length-diameter ratio of 2, and a volume of 1270 cm 3 ; 

C - a similar Centaur liquid-hydrogen-tank model that included internal 
slosh baffles at the 30- and 77-percent liquid fill levels; and D - a 
spherical tank with a volume of 524 cm 3 . 

The piston pump was located in the convex hemispherical end of each 
tank. This location provided simulation of the liquid flow patterns 
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TABLE I. - SUMMARY OF TEST CONDITIONS AND ANALYTICAL RESULTS 
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associated with a pump in the forward end of each tank. The forward end of 
the tank was chosen as the pump location because under Earth-orbiting con- 
ditions liquid would be positioned over the pump as a result of the acceler- 
ation caused by atmospheric drag. 

For this experimental study, tank Bond numbers (acceleration to surface 
tension force ratio) of 2, 4, 8, and 12 were chosen to cover the range of 
accelerations expected on orbit transfer vehicle propellant tanks. These 
Bond numbers correspond to accelerations on the experiment package of 
0.0023, 0.0046, 0.0092, and 0.014 g's, which were achieved by the use of a 
cold gas thruster. Other experimental variables were initial experiment 
tank liquid filling, axial liquid-jet outlet position, and liquid-jet 
volumetric flow rate. 

The lengths of the axial jet outlet tubes, 1 and 8 cm, were measured 
from the hypothetical bottom of each experiment tank. For each test the 
average piston pump velocity was determined from the position-versus-time 
trace obtained from the strip-chart recorder. The average piston pump 
velocity was multiplied by the piston cross-sectional area to obtain the 
liquid-jet volumetric flow rate. 


EXPERIMENTAL RESULTS 

The experimental conditions for the 14 normal-gravity tests, 36 reduced- 
gravity tests, and 41 zero-gravity tests (<10 - ^ g's), together with the 
characteristics of the axial liquid jet for each test, are presented in 
table I. The data from all of the reduced-gravity tests and those from the 
zero-gravity tests with initial experiment-tank liquid fillings of approxi- 
mately 40 and 60 percent (by volume) have not been previously published; the 
remaining data are taken from references 4 and 5 for comparison. 

Four distinct liquid flow patterns, previously observed and discussed in 
reference 5, are shown in figure 3 and identified by the following Roman 
numerals in table I: I - dissipation of the liquid jet in the bulk-liquid 

region (fig. 3(a)), II - geyser formation (fig. 3(b)), III - collection of 
jet liquid in the aft end of the tank (opposite the jet outlet) (fig. 3(c)), 
IV - liquid circulation over the aft end of the tank and down the tank wall 
(fig. 3(d)). 

Many of the tests exhibited results that combined characteristics of two 
flow patterns. Experimental conditions that produced some bulk-liquid mix- 
ing together with the formation of fairly prominent geysers are identified 
by the Roman numerals I / 1 1 . Tests identified as having flow patterns 1 1 / 1 1 1 
or II/IV indicate the author's uncertainty due to the limited test time. A 
few tests that resulted in liquid flow patterns that do not fall into any of 
the previously discussed categories are identifed by footnotes. 

The liquid height above the liquid-jet outlet was measured just prior to 
the movement of the piston pump. The maximum height of each geyser was also 
measured and divided by the tank radius. The time required, after the pis- 
ton pump began to move, for complete bulk-liquid mixing was either measured 
or extrapolated from the data for each test. The measured bulk-liquid mix- 
ing times are the times required for the bulk liquid to become completely 
dyed by the incoming liquid jet. The shortcomings associated with the use 
of a dyed liquid jet in an isothermal system to simulate the mixing of 
noni sothermal liquids are examined in the section Concluding Remarks. The 
extrapolation technique involved measuring the position of the demarcation 
line between the dyed and undyed liquids in the bulk-liquid region as a 
function of time, from these measurements the velocity of the demarcation 
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(a-1) Normal gravity, test 16. (a-2) Weightlessness, test 17 
ia) Dissipation of liquid jet in bulk-liquid region. 

(b) Geyser formation, test 24. 

(c) Liquid collection, test 18. 

Id) Liquid circulation, test 42. 

Figure 3. - Observed axial liquid-jet flow patterns. 




line between the dyed and undyed liquids was determined and multiplied by 
the tank cross-sectional area to establish the bulk-liquid mixing rate. The 
volume of unmixed liquid bulk at the end of each test was then divided by 
this experimentally determined mixing rate, and the resulting time increment 
was added to the test time to obtain an estimate of the time required for 
complete bulk mixing. Tests for which bulk-liquid mixing times are not 
presented fall into three categories: 

(1) Low liquid-jet flow rates coupled with only 3 sec of test time yield 
immeasurable rates of bulk-liquid mixing (flow patterns I and IV). 

(2) Mixing of the liquid jet and the bulk liquid is negligible (flow 
patterns 1 1 and III). 

(3) Bulk-liquid mixing phenomena are affected by the presence of slosh 
baffles. 


ANALYSIS OF EXPERIMENTAL RESULTS 

Dimensional Analysis of Liquid-Jet Flow Patterns 

Studies of liquid inflow to model tanks under reduced-gravity conditions 
(refs. 6 and 7) provided the basis for the author's analysis of the experi- 
mental data presented herein. The height of the liquid geyser that forms 
during inflow was selected as the most prominent physical characteristic of 
the observed fluid dynamic phenomena. The geyser height is also a suitable 
physical characteri Stic to describe three of the four flow patterns that 
were observed during the propellant mixing studies: flow pattern I is dis- 

tinguished by negligible geyser heights, flow pattern II is described 
directly by geyser height, and flow pattern III is characterized by geysers 
that just touch the aft end of the tank opposite the jet outlet. 

The inertia, acceleration, and surface tension forces are the primary 
factors that determine geyser formation. Examination of the formulas that 
describe these forces leads to the selection of the liquid-jet velocity, 
density, and surface tension, the acceleration, the geyser height, and 
several characteristic dimensions as the important parameters in this study. 

The Buckingham theorem was employed to establish three dimensionless 
groups that can be used to correlate the experimental data: the Bond Bo 

and Weber We numbers and a nondimensional geyser height G. The Bond 
number is the ratio of acceleration to surface tension forces, and the Weber 
number is the ratio of inertial to surface tension forces. To avoid 
confusion, the reader should recognize the difference between the liquid jet 
Bond number and the tank Bond number. The liquid-jet Bond number, based on 
the jet radius Rj, is used to characterize the interaction between the 
liquid jet and the liquid-vapor interface. The tank Bond number, based on 
the tank radius R^, is used to characterize the bulk-liquid behavior. 

The three dimensionless groups are 


G = 


h/R. 
9 t 


( 1 ) 


10 


Bo = apR ^ / a 

J 


( 2 ) 



( 3 ) 


We = pvV/aD . 

0 0 J 

With one exception the parameters contained in these three dimensionless 
groups are either fluid properties, easily measured dimensions, or terms 
that can be readily calculated by using the continuity equation. Establish- 
ing the functional relationship between the experimental variables and the 
radius of the liquid jet proved to be a significant portion of the 
analytical effort. 

Axial liquid-jet physical characteristics . - The velocity with which the 
liquid jet exits the simuTated pump affects not only the momentum of the jet 
and thus the resulting liquid flow pattern, but also the physical appearance 
or characteristics of the liquid jet (ref. 8). At very low velocities the 
liquid jet is laminar, with almost undetectable spreading of the jet as it 
passes through the bulk liquid (fig. 3 ( a— 1 ) j . As the liquid jet velocity is 
increased, a transition to completely turbulent flow occurs, and the jet 
spreads into the bulk region at an angle of up to 12° on each side of the 
jet centerline (fig. 3(b)). 

Because of the liquid-jet spreading phenomenon, the position of the 
outlet for an axial liquid jet relative to the bulk liquid-vapor interface 
greatly affects the resulting jet characteri sties. According to the 
classical theory for the behavior of axial jets flowing through a fluid 
(ref. 9), a submerged jet maintains nearly constant momentum as it 
progresses from the outlet. As the jet flows through the bulk liquid, its 
mass flow rate increases because of bulk-liquid entrainment, and its average 
velocity decreases because of radial momentum transport by viscous effects. 
Consequently, for axial jets with outlets positioned in the bulk liquid, 
increasing the liquid depth over the axial liquid-jet outlet will increase 
the size and mass flow rate of the liquid jet and reduce its average 
velocity at the liquid-vapor interface. Comparison of figures 3(b) and (c) 
provides dramatic evidence of this fluid dynamic phenomenon. 

The jet outlet position and tank filling are not the only parameters 
that determine the depth of bulk liquid that the liquid jet must penetrate. 
Under weightless conditions the liquid-vapor interface becomes hemis- 
pherical, and thus liquid depth over the jet outlet is less than with the 
same tank filling under normal gravity. Two exceptions to t his occur in the 
baffled Centaur model tank. For tank liquid fillings of approximately 30 
and 70 percent the lower and upper slosh baffles, respectively, restrict the 
motion of the liquid-vapor interface so that the highly curved zero-gravity 
configuration typical of unbaffled tanks does not occur. Consequently, the 
liquid depth over the jet outlet was greater for these test conditions than 
for the corresponding tests in the unbaffled tanxs. 

Correlation of liquid- jet characteristics . - Reference 6, which presents 
a summary of the pertinent literature, suggests that, at liquid-jet Reynolds 
numbers greater than approximately 1500, the jet should be completely tur- 
bulent. The following form of the Reynolds number was used for this study: 


Re 


pv D /u 
oo 


(4) 


The calculated value of the Reynolds number for each experimental test can 
be found in table I. 
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A completely turbulent jet should be expected to spread to about 6 
to 8° on each side of the get centerline for a distance of 12.4 jet outlet 
radii and then to about 10 to 12° for greater distances, from these angles 
of spread, the following approximate expressions for the radius of the 
liquid jet can be derived (ref. 6), where h b is the distance between the 
jet outlet and the liquid-vapor interface: 


R. = R + 0.12 h. h. < 12.4 R (5) 

j,c o b b - o 


R . = 0.11 R + 0.19 h, h. > 12.4 R (6) 

j,c o b b o 

For laminar liquid jets the angle of spread is only 2 to 3 on either side 
of the jet centerline, for these conditions of jet spreading the approxi- 
mate expression for the radius of the liquid jet as a function of the liquid 
height above the liquid-jet outlet is 



= R + 0.04 h. 
o b 


(7) 


The angle of spread of the liquid jet, based on measurement of the 
penetration of the dyed fluid into the bulk-liquid region, was determined 
for each experimental test. Comparison of the experimentally obtained 
measurements of the liquid-jet dimensions with the previously presented 
theory led to the following conclusions: 

(1) For experimental conditions that yielded Reynolds numbers in excess 

of 1250 and liquid heights above the liquid-jet outlet greater than 2.5 cm 
(12.4 R 0 ) , agreement between equation (6) and the experimental results was 
good. 

(2) The data obtained at a liquid-jet Reynolds number of 630 and liquid 
heights greater than 2.5 cm agreed almost exactly with equation (7). 

(3) At liquid-jet Reynolds numbers of approximately 450 no spreading of 
the jet could be measured. 

(4) When the liquid-jet outlet was less than 2.5 cm below the liquid- 
vapor interface, the resulting jet had the same radius as the jet outlet 
tube. This experimental observation conflicts with the predictions of 
equation (5); the reason for the difference could not be determined. 

(5) For experimental conditions that yielded Reynolds numbers of 
approximately 900 and liquid heights greater than 2.5 cm, the liquid-jet 
characteristics were intermediate between the laminar and completely tur- 
bulent results. The following equation was found to fit the experimental 
data for these conditions: 


R . = 0.10 R + 0.13 h 

J,P o 


h, > 12.4 R 
b o 


(8) 


The preceding criteria were used to determine the liquid-jet radii at the 
liquid-vapor interface for each test. Then, by using equations (2) and (3), 
the liquid-jet Weber and Bond numbers were calculated and recorded in table I. 



In addition, for those tests that resulted in flow pattern II, III. or 
the width of the geyser was measured just above the liquid-vapor interface. 

Examination of these data led to the following conclusions: 

(1) When the liquid-jet outlet was less than 2.5 cm below the liquid 
vapor interface, no dissipation of the liquid jet in the bulk liquid was 
apparent, and the resulting geyser had the same radius as the jet out e 

tUbe '(2) For liquid heights above the liquid-jet outlet greater than 2.5 cm, 
the radius of the resulting geyser was 10 to 20 percent smaller than the 
radius of the liquid jet at the liquid-vapor interface. This observation 
can be explained by realizing that the geyser is formed from the higher 
velocity central core of the liquid jet while the lower velocity outer 
extremities of the jet are dissipated in the bulk liquid. 

Correlation of lio uid-jet dimensionless groups . - To establish the 
functional relationship between the 'three dimensionless groups established 
by the Buckingham theorem, all of the partially turbulent (Reynolds number, 
aDDrox 900) zero-gravity cylindrical tank data were examined. This initial 
screen i ng , which Included only zero Bond number data, led to the conclusum 
that a linear functional relationship exists between the nondimensional 
geyser height and the liquid-jet Weber number: 

G = f(X + V We) ^ 

The logical assumptions that the geyser height should be inversely pro- 
portional to the magnitude of the acceleration environment and that t e 
correlating expression must reduce to equation (9) under zero-gravity 
conditions yield 


G = f 


X + Y We 
1 + Z Bo 


Use of a least-squares deviation curve-fitting technique for the 
partially turbulent zero-gravity experimental data established values for 
the constants X and Y of -0.3 and 1.0, respectively. The negative value for 
X reflects the fact that a finite jet velocity and corresponding Weber num- 
ber are required before the liquid jet is not totally dissipated in e 
liquid bulk and a geyser forms. Further curve fitting of the reduced- and 
normal-gravity data established a value of 0.6 for the constant Z. Thus 
flow characterization parameter Fn was established that . can be used o 
describe the flow patterns resulting from partially turbulent (h b > 2.5 cm) 
submerged liquid jets with Reynolds numbers of approximately yuu: 


_ We - 0.3 
p " 1 + 0.6 Bo 


Figure 4(a) shows the excellent correlation that was achieved between 

the experimentally determined nondimensional geyser heights and the cal- 
culated flow characterization parameters for the partially turbulent liquid 

jets. 
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A correlating equation was similarly derived for completely turbulent 
(Re > 1250) submerged liquid jets: 


F _ 1.6 We - 0.5 
c 1 + 0.6 Bo 


( 12 ) 


Figure 4(b) is a plot of the experimental ly determined nondimensional 
geyser height versus the flow characterization parameter for the completely 
turbulent liquid jets. The data scatter exhibited by this figure, as com- 
pared with figure 4(a), should be expected based on the wide range of flow 
rates that are included in the turbulent liquid-jet data. Also, the start- 
ing points for this analysis, reference 6 and equations (5) and (6), were 
established based on steady-state experimental conditions for a gas jet 
flowing into a gaseous environment. Although gas-gas and liquid-liquid 
systems should exhibit similar characteristics (ref. 9), the short test 
times associated with this experimental study suggest that transient effects 
could also contribute to the observed data scatter. 

Values of the flow characterization parameter for both the partially and 
completely turbulent liquid jets are recorded in table I. Because of the 
limited amount of laminar liquid-jet data no attempt was made to correlate 
the experimental results for these test conditions. However, a comparison 
of the Weber numbers and the nondimensional geyser heights for the laminar 
liquid-jet tests indicates that the same general functional relationship 
exists. For the test conditions where the liquid jets exited near the 
liquid-vapor interface (h b < 2.5 cm), only flow patterns III and IV 
resulted under zero or reduced-gravity conditions. Consequently, no cor- 
relation of the experimental results was possible since the qeysers or iets 
are constrained by the test container and their maximum height cannot be 
determined. 

Flow characterization parameters were also calculated by using equa- 
uons (11) and (12) for the tests conducted in spherical containers 
(table I). Although the same general trends exist, good agreement between 



F . (We -013) 
P (0.6 Bo + 11 



F _ (1. 6 We - ft 5) 
c (ft 6 Bo +• 1) 


(a) Jet Reynolds number, Rej = 900. (b) Jet Reynolds number, Rej > 1250. 

Figure 4. - Nondimensional geyser height as a function of flow characterization parameter. Cylindrical tank data only. 


the experimentally determined nondimensional geyser height and the cor- 
responding values of the flow characterization parameter did not result; a 
possible explanation for this variance is presented in the section 
Discussion of Results. 


Bulk-Liquid Mixing Correlation 


References 10 and 11, which present an extensive review and summar- 
ization of fluid mixing theory, were the sources of several liquid-mixing 
correlating parameters. Values of these correlating parameters were 
calculated for the test results reported herein for which bulk-liquid mixing 
times and flow characterization parameters had been determined. A mixing 


correlation parameter eQh^D^y originally developed by Okita and 

Oyama (ref. 12) was found to be a function of the flow characterization 
parameters over a wide range of experimental variables. Since the work of 
Okita and Oyama dealt only with completely filled tanks, the bulk-liquid 
volume Vb was substituted for the tank volume Vt to account for the 
variation in the liquid filling that was one of the variables in this 
exper i ment a 1 program. In addition, the characteristic dimensions contained 

in the reference 12 mixing parameter were either constants (D^ 2 and Dj) or 

varied only slightly (hj^ 2 ) in this experimental program. Consequently, the 


Okita and Oyama mixing parameter was reduced to a much simpler form - the 
liquid-jet volumetric flow rate divided by the bulk-liquid mixing rate: 


In this expression the liquid-jet volumetric flow rate Q is an 
independent variable. The bulk-liquid volume V^ is dependent on Doth 
the initial tank filling and the liquid-jet characteristics. The bulk- 
liquid mixing time 0 is a dependent variable whose values were 
determined experimentally. To make the experimental results more useful, an 
analytical technique was developed that provides a means of determining 
approximate values for the bulk-liquid volumes. 

The volume of liquid initially in the experiment tanks is reduced by the 
action of the liquid jets through the formation of geysers (flow patterns 
II, III, and IV) and the presence of liquid on the tank walls in the end of 
the tank opposite the mixer unit (flow patterns III and IV). Since no 
bulk-liquid mixing was observed for test conditions that yielded flow 
pattern III, the analysis was restricted to examination of flow patterns II 

and IV. . 

For test conditions that produce only geysers the bulk-liquid volume can 
be determined by subtracting the volume of liquid in the geyser from the 
initial volume of liquid in the experiment tanks; the difference between the 
liquid-jet volumetric flow rate and the rate of liquid withdrawal at the 
base of the simulated mixer unit is small and was neglected. The volume of 
liquid in the geysers can be approximated by using the previously presented 
analytical techniques. On the basis of the liquid jet Reynolds number the 
appropriate equations for the liquid-jet radius at the liquid-vapor inter- 
face, equation (6) or (8), and for the nondimensional geyser height, 


equation (11) or (12), can be employed to calculate the bulk-liquid volume 
(flow pattern II): 


V b = V t (% filling) - ifRjFR t 


(14) 


For test conditions that produce flow pattern IV both the volume of 
liquid in the geyser and the volume on the tank walls must be subtracted 
from the initial volume of liquid in the experiment tank to establish the 
bulk-liquid volume. The geyser radius is once again approximated by usinq 
either equation (b) or (8), and the geyser height is now the distance from 
the liquid-vapor interface to the end of the tank opposite the mixer unit- 


h g - h t ~ h b " h o (15) 


Analytical techniques originally developed by Prandtl and subsequently 
experimentally verified (ref. 13) indicate that liquid jets which impinge on 
a solid surface can be characterized as having uniform and continuous veloc- 
ity profiles. The only exceptions to this general observation occur in the 
thin boundary layer of fluid adjacent to the tank wall and in the immediate 
vicinity of the point of impingement of the liquid jet. Consequently for 
the engineering solution required by this analysis, it is possible to’ 
neglect the boundary layer and jet impingement point effects and to assume 
that the average velocity of each geyser is equal to the resulting averaqe 
velocity of the liquid flowing on the tank walls. This assumption is 
reasonably applied only to the hemispherical ly ended tanks (tank shapes 
A and D). However, this restriction is of no consequence since no bulk- 
liquid mixing times were established for experimental conditions that 
yielded flow pattern IV in the Centaur model tanks (tank shapes B and C). 
Because continuity must be preserved as the liquid flows on the tank wall 
^oss-sectional area of this liquid is constant and only the thickness 
or the liquid changes to accommodate the geometry of the tank. For tank 
shapes A and 0 a simple geometric analysis yields the following expression 
for the bulk-liquid volume (flow pattern IV): 


V b = V t^ flll ' in 9) - *R 2 

J 


The values of the mixing parameter M (eq. (13)) calculated by using 
bulk-liquid volumes determined from equation (14) (flow pattern II) or 
equation (16) (flow pattern IV) can be found in table I. 

Figure 5 presents a plot of the mixing parameter as a function of the 
flow character i zati on parameters for the data obtained in the cylindrical 
tanks with a jet outlet tube 1 cm long. The data obtained with partially 
turbulent liquid jets generally exhibited more scatter than data for com- 
pletely turbulent liquid jets. This observation can be explained by 
realizing that the partially turbulent liquid jets result from relatively 
low liquid- jet flow rates. Because of the limited test time available, low 
liquid-jet flow rates yield the least observed quantity of bulk-liquid 
mixing and thus the greatest extrapolation of the experimental data. 


2h + R, 


(t- 1 ) - R j 


(16) 
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Figure 5. - Mixing parameter as a function of flow characterization parameter. Cylindrical tanks with 
1-cm jet outlet tubes. 


For experimental conditions that yield liquid flow patterns I and IV, 
the mixing parameter increases steadily as the flow characterization paran>- 
eter increases. Although these flow patterns are strikingly different, the 
resulting bulk-liquid mixing phenomena are both observational ly and analyti- 
cally similar. The following functional relationship exists between the 
mixing and flow characterization parameters for liquid flow patterns I 
and IV: 


M = 0.09 + 0.01 In F 


(17) 


In the transition region from flow pattern I to flow pattern II, there 
is a lot of scatter in the values of the mixing parameter. This data scat- 
ter could be anticipated from observing the flow phenomena, which exhibit a 
great deal of irregular flow and a pulsing motion of the geysers that form. 
This pulsing motion aids bulk-liquid mixing so that the mixing parameter 
decreases in value as the flow characterization parameter increases. As the 
flow characterization parameter increases to the point where large geysers 
result (flow pattern II), the bulk-liquid mixing times are longer because 
mixing is negligible. For experimental conditions that yield flow pattern 
III, no bulk-liquid mixing occurs and the value of the mixing parameter 
approaches infinity. 

Values of the mixing parameter were also calculated for the experimental 
tests conducted in the spherical tanks and for the one test (no. 63) that 
was conducted in a cylindrical tank with an 8-cm-long liquid-jet outlet 
tube. The same general trends in the functional relationship between the 
mixing parameter and the flow characterization parameter were observed for 
these data. However, the mixing parameters were generally lower at the same 
flow conditions. Because of the uncertainty in the validity of applying the 
f low-characterization-parameter calculation procedure for the 1-cm-long- 
1 iquid-jet-outlet-tube cylindrical tank to the spherical tank and test 63 
data, these mixing data were not plotted in figure 5. 


DISCUSSION OF RESULTS 


A few data points included in table I exhibited liquid-jet flow 
characteristics and bulk-liquid mixing phenomena that do not fit into the 
general categories introduced in the previous sections. 

Test 56 was the only data point obtained for both an 8-cm-long outlet 
tube and appreciable liquid depth at a jet Reynolds number of approximately 
900. Observation of the penetration of the liquid-jet dye into the bulk 
liquid revealed a thin laminar jet similar to the liquid jets that resulted 
at lower Reynolds numbers with the 1-cm-long outlet tube. It is the 
author's opinion that the longer jet outlet tube yields a more completely 
developed liquid-jet velocity profile and that the transition to partially 
turbulent flow thus does not occur. Consequently, the Weber and Bond num- 
bers for test 56 were calculated by using equation (7) to determine the 
radius of the laminar liquid jet. Test 63, which had both an 8-cm-long 
outlet tube and appreciable liquid depth at a jet Reynolds number of 1850, 
was observed to have a completely turbulent liquid jet. 

Tests 38, 45, and 69 were all conducted under normal-gravity conditions 
with liquid depths over the liquid-jet outlet tube of approximately 2.5 cm 
and jet Reynolds numbers that ranged from 900 to 5150. For these test con- 
ditions the geysers that formed had the appearance of a "fountain" with no 
distinct liquid-vapor interface (footnote a in table I), and the incoming 
liquid fell back into the bulk-liquid region, causing centralized mixinq 
only. 

At the highest liquid-jet flow rates tested, Reynolds numbers of approx- 
imately 7500, undesirable liquid flow characteristics resulted under both 
zero- and normal-gravity conditions (tests 46, 48, and 49). Under zero- 
gravity conditions the liquid flow down the tank walls was entrained by the 
liquid jet so that most of the dyed liquid recirculated in the vapor region 
of the tank and very little bulk-liquid mixing resulted (footnote b in 
table I). For all three tests these high liquid-jet flow rates caused 
bubbles to be entrained in the bulk liquid (footnote c in table I) and thus 
defeated the objective of providing bulk liquid of uniform temperature and 
density. 

Spherical tank tests 81 and 83, conducted under zero-gravity conditions 
with liquid fillings of approximately 70 percent, exhibited bulk-liquid col- 
lection on the sides of the container after the start of liquid-jet flow. 

The jet liquid passed through the central portion of the container and 
flowed over the bulk liquid to the bottom of the tank (footnote d in 
table I); no bulk-liquid mixing resulted. At experimental conditions of 
zero gravity and approximately 90 percent liquid filling of the spherical 
tank (tests 87 and 89), the start of liquid-jet flow caused the bulk liquid 
to collect at the end of the tank opposite the jet exit and the vapor bubble 
to be drawn into the simulated mixer inlet (footnote e in table I). 

The author anticipated that the application of a small acceleration 
would be sufficient to position the bulk liquid positively and thus to avoid 
the just discussed problems associated with relatively high liquid fillings 
in spherical tanks. In general, the stabilizing influence of a small 
acceleration did in fact lead to desirable liquid flow patterns and uniform 
mixing of the liquid bulk. However, test 90 proved to be the exception that 
leads to the conclusion that axial liquid-jet mixing of propellants in 
spherical tanks will still be difficult even under reduced-gravity con- 
ditions. For this test, flow pattern IV was observed at a liquid filling of 
approximately 90 percent, but the bulk-liquid mixing was restricted to the 
central portion of the container (fig. 6). 
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Figure 6. - Low-gravity bulk-liquid mixing; liquid filling, 
90 percent; test 90. 


Application of the cylindrical tank flow-characterization-parameter 
calculation procedure to the spherical tank data yielded less than satis- 
factory results. The calculated flow characterization parameters were 
consistently higher than the experimentally determined nondimensional geyser 
heights. Under zero- and reduced-gravity conditions the liquid-vapor inter- 
face has a smaller radius of curvature in spherical tanks than in cylindri- 
cal tanks because of the curvature of the tank wall. The author suggests 
that, because of this more highly curved liquid-vapor interface, the radius 
of the spherical tank is probably not an appropriate characteristic dimen- 
sion to be used to nondimensional ize the geyser height. However, since the 
amount of spherical tank data was limited, no attempt was made to establish 
a more appropriate characteristic dimension or flow characterization param- 
eters specifically for spherical tanks. 


APPLICATION OF MODELING TECHNIQUE 

To demonstrate the use of the correlations developed in the section 
Analysis of Experimental Results, a typical liquid-hydrogen tank and 
internal thermodynamic vent system designed for space operation have been 
analyzed. The selected propellant tank (ref. 14) and internal thermodynamic 
vent system (ref. 3) are for a shuttle-compatible space tug that is very 
similar to current orbit transfer vehicle concepts. Table II is a summary 
of the parameters used for this analysis and the resulting recommendations. 

It was assumed that the liquid-hydrogen tank is covered with multilayer 
insulation (30 layers) so that the resulting in-space heating rate of 
approximately 88 W (300 Btu/hr) would match the cooling capability of the 
reference 3 thermodynamic vent system. At this design point the mixer unit 
provides a liquid-hydrogen flow rate of 2300 cnr/sec (1320 lb/hr) through 
the 4.4-cm (1.75-in.) diameter outlet tube. The corresponding overboard 
vent rate is 0.8 kg/hr (1.7 lb/hr). The electric power required for the 
mixer at this design point is approximately 7 W, with a resulting increase 
in the rate of heat addition to the tank contents of 2.3 W (8 Btu/hr). 
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TABLE II. - APPLICATION OF MODELING 
TECHNIQUE PARA INTERS AND RESULTS 
Liquid-hydrogen tank: 

Diameter, m(ft) 4.3 (14) 

Length, ml ft) 6.1 (20) 

Volume, m J (ft J ) 71 (2500) 

Heating rate, W(Btu/hr) 88 (300) 

Reference 3 mixer unit characteristics: 

Flowrate, cr^/secj lb/hr) 2300 (1320) 

Outlet tube diameter, cm(in.j 4.4 (1.75) 

Cooling rate, W(Btu/hr) 88 (300) 

Overboard vent rate, kg/hr( 1 b/hr ) 0.8 (1.7) 

Recommended mixer unit and resulting 
liquid-jet characteristics: 

Flowrate, cm^/sec(lb/hr) 3450 (1980) 

Outlet tube diameter, cm(in.j 17 (6.7) 

Cooling rate, W(Btu/hr) 132 (450) 

Overboard vent rate, kq/hr( lb/hr) 1.2 (2.6) 

Jet velocity, cin/sec(ft/sec) 15 (3.2) 

Jet Weber number 3 3.3 to 12.1 

Jet Bond number 3 0.85 to 0 

Flow characterization parameter 3 .... 3.2 to 18.9 


Tank fillings range from 90 to A 20 vol* and acceleration 
environment ranges from 3xl0" b g's to 0. 


A typical multiburn orbit-raising mission was selected for the orbit 
transfer vehicle. Consequently the liquid-hydrogen propellant tank 
conditions varied from, initially, a nearly full tank in low Earth orbit 
withan acceleration level (caused by atmospheric drag) of approximately 
3xl0~ b g to, just before the final engine burn, a tank approximately 20 
percent full and an acceleration level of essentially zero. 

Because of the low viscosity of liquid hydrogen the liquid-jet Reynolds 
number is very high and equation (6) can be used to establish the radius of 
the completely turbulent liquid jet at the liquid-vapor interface. Weber 
numbers were calculated for a range of orbit transfer vehicle hydrogen-tank 
liquid fillings by using equation (3). At the mixer unit design-point flow 
rate of 2300 cnr/sec the liquid-jet Weber number exceeds 60 at low tank 
fillings. On the basis of model tank experimental data a submerged liquid- 
jet Weber number in excess of 30 indicates that turbulent, recirculating 
liquid flow plus the formation of bubbles will probably result - a highly 
undesirable situation. 

Increasing the diameter of the liquid-jet outlet tube while maintaining 
the same mixer flow rate and thus the same heat-exchanger cooling capacity 
will reduce the liquid-jet Weber numbers to within acceptable limits. A 
liquid-jet outlet tube diameter of 17 cm (6.7 in.) was chosen to provide 
dimensional similitude with the model experiment tanks (tank diameter to jet 
outlet diameter ratio, 25). The liquid-jet Weber and Bond numbers were 
calculated by using equations (2) and (3) with a resulting maximum value for 
the Weber number of 5. The corresponding values of the flow character- 
ization parameter varied from 1.2 to 7.5. Consequently, as the tank filling 
decreases during the mission, a transition from flow pattern II to flow pat- 
tern IV would occur - a highly undesirable situation. 

To achieve dissipation of the liquid jet in the bulk-liquid region (flow 
pattern I), the flow characterization parameter should be less than 0.5. 
Obviously, for the orbit transfer vehicle chosen as an example, achieving 
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flow pattern I would be impossible without either using a very large- 
diameter jet outlet tube or reducing the liquid- jet flow rate and thus the 
cooling capacity of the thermodynamic vent system. Of more concern is the 
avoidance of flow patterns II and III and the resulting lack of bulk-liquid 
mixing. The author suggests that, in addition to using a 17-cm-diameter 
liquid-jet outlet tube, the jet flow rate should also be increased by 
approximately 50 percent. This higher liquid-jet flow rate will yield jet 
Weber numbers and flow characterization parameters that will guarantee 
complete liquid circulation (flow pattern IV) without the possibility of 
undesirable liquid flow conditions. In addition, the higher flow rate will 
provide additional thermodynamic vent system cooling capacity to allow for 
any uncertainty in the rate of heat input to the liquid-hydrogen tank. 

Further examination of orbit transfer vehicle requirements reveals that 
the heating load on the liquid-hydrogen tank while it is still in the 
shuttle cargo bay may be as much as five times higher than the in-space 
heating load (ref. 15). This higher heat transfer rate results primarily 
from the time required for the multilayer insulation to become evacuated 
following exposure to the vacuum of space. If an internal thermodynamic 
vent system must provide 1 iquid-hydrogen-tank pressure control while the 
orbit transfer vehicle is in the shuttle cargo bay, a much larger system 
than that reported on in reference 3 will be required. To accommodate the 
complete range of propellant tank conditions (heating rate, liquid filling, 
and acceleration environment), it may be necessary to design the thermo- 
dynamic vent system for both intermittent operation and variable mixer unit 
flow rate capability. 


CONCLUDING REMARKS 
Fluid Dynamic Analysis 

The analytical technique of employing dimensional analysis to provide 
fluid dynamic criteria for the in-space operation of liquid-fueled propul- 
sive stages is well established (ref. 16). Specifically, the use of the 
Bond and Weber number nondimensional parameters to establish the similarity 
between ground-based experimental programs conducted in drop towers using 
scale-model tanks and the in-space operation of full-scale vehicles has been 
demonstrated (ref. 17). This analytical technique has been applied to the 
problems of establishing the criteria for liquid-vapor interface stability, 
liquid reorientation, propellant tank draining, and liquid inflow to propel- 
lant tanks under reduced-gravity conditions. 

The fluid dynamic analysis presented herein is a logical extension of 
the dimensional analysis technique to the problem of mixing propellants in 
space. For the cylindrical tank data an excellent correlation resulted 
between the experimentally determined nondimensional geyser heights and the 
newly defined flow characterization parameter (fig. 4). That this corre- 
lation was good over a wide range of liquid-jet flow conditions and gravi- 
tational environments lends credence to the use of this technique to predict 
axial liquid-jet flow patterns in full-scale propellant tanks. 


Fluid Mixing Analysis 

Liquid-jet mixing parameters have been shown to be a function of the jet 
Reynolds number under normal gravity (refs. 10 and 11). A previous attempt 


(ref. 18) to correlate reduced- and normal-gravity liquid-jet mixing 
phenomena also used the jet Reynolds number to characterize the liquid jet. 
This approach led to vastly different values for the mixing parameters when 
experimental data with the same jet Reynolds number but different acceler- 
ation environments were compared. 

The most significant contribution to liquid mixing theory reported 
herein is the establishment of flow characteri zation parameters that are 
functionally related to a simplified version of the Okita and Oyama mixing 
parameter (ref. 12) over a wide range of experimental conditions. Of par- 
ticular importance is the fact that the range of acceleration environments 
explored covered several orders of magnitude. 

The use of dyed ethanol to simulate the mixing of cryogenic propellants 
does, however, leave two unresolved issues: (1) how the mixing of dyed and 

undyed liquids is related to the mixing of liquids of different temperature; 
and (2) how the differences in fluid properties between ethanol and typical 
cryogenic propellants influence the mixing parameter. The correlation be- 
tween the flow characterization parameters and the simplified Okita and 
Oyama mixing parameter that has been established suggests that future 
experimental data obtained under normal-gravity conditions can be used to 
resolve these issues and to predict the in-space performance of axial 
liquid-jet mixing devices. 


SUMMARY OF RESULTS 

An experimental program was conducted in the Lewis Research Center's 
zero-gravity facility to examine the liquid flow patterns that result from 
the axial jet mixing of ethanol in 10-cm-diameter containers under zero-, 
reduced-, and normal-gravity conditions. Four tank configurations were used 
for the study: a spherical container; a convex hemispherically ended cylin- 

drical container; and two Centaur liquid-hydrogen-tank models, one with 
slosh baffles and one without. 

How the gravitational environment, the tank configuration, the position 
of the liquid-jet outlet tube relative to the liquid-vapor interface, and 
the liquid-jet volumetric flow rate affected the liquid flow patterns was 
examined. Four distinct liquid flow patterns were observed: 

(1) Pattern I - dissipation of the liquid jet in the bulk-liquid region 

(2) Pattern II - geyser formation 

(3) Pattern III - collection of the jet liquid in the aft end of the tank 
(opposite jet outlet) 

(4) Pattern IV - liquid circulation over the aft end of the tank and 
down the tank wall 

Nondimensional parameters (Reynolds, Weber, and Bond numbers) that are 
typically applied to the analysis of fluid dynamic phenomena were used to 
characterize the liquid jets and the resulting liquid flow patterns. The 
Reynolds number provides differentiation between the laminar, transition, 
and turbulent liquid-jet flow regimes. Tne experimentally observed geyser 
height, divided by the tank radius, was chosen to characterize the observed 
flow patterns. This nondimensional geyser height was correlated with a new 
nondimensional grouping of the Weber and Bond numbers (flow characterization 
parameters) in both the transition and turbulent liquid-jet flow regimes. 

The liquid-bulk mixing phenomena, expressed as the ratio of liquid-jet 
volumetric flow rate to liquid-bulk mixing rate, for liquid flow patterns I 
and IV were shown to be a function of the flow characterization parameters 
over a wide range of experimental variables. Although a great deal of 
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scatter exists in the data, this mixing parameter also adds insight to the 
liquid-bulk mixing phenomenon for liquid flow pattern II, No liquid-bulk 
mixing occurs for liquid flow pattern III. 

Use of the flow characterization parameters to determine mixer require- 
ments for a full-scale orbit transfer vehicle (OTV) hydrogen tank yields the 
following conclusions: 

1. A mixer designed to achieve dissipation of the liquid jet in the 
bulk-liquid region (flow pattern I) will require low liquid-jet flow rates 
and large heat exchangers. 

2. The recommended mixer design is based on achieving complete liquid 
circulation (flow pattern IV) during the OTV orbit-raising mission. 

3. The anticipated variation in tank heating rates for cryogenic tankage 
transported to orbit by the shuttle may require multiple flow rate capa- 
bility and/or periodic operation of the mixer unit. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, July 2, 1982 
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